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Abstract: The importance of direct capture for (n,7)-reactions on intermediate- 
and heavy-mass target nuclei occuring in the s- and r-process is investigated. It 
is shown that the direct mechanism is non-negligible for magic and neutron rich 
target nuclei. For some double magic and neutron rich nuclei in the r-process direct 
capture is even the dominant reaction mechanism. 
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1 Introduction 



The reaction mechanism of radiative capture can be classified by two extremes: the 
reaction proceeds as a multi-step process (compound-nucleus reaction: CN) or in 
a single step (Direct Capture: DC). In the last years it was realized that for light 
nuclei DC is often the dominating reaction mechanism in astrophysically relevant 
nuclear processes. The DC-method together with the folding procedure has been 
applied sucessfully by our group to radiative capture in primordial nucleosynthe- 
sis ( 2 H(a, 7 ) 6 Li 0, 3 H(a, 7 ) 7 Li |], 7 Li(n, 7 ) 8 Li §), the pp-chain ( 3 He(a, 7 ) 7 Be §, 
7 Be(p, 7 ) 8 B §), CNO-cycle ( 13 N(p, 7 ) 14 §) and helium burning ( 8 Be(a, 7 ) 12 C [§). 

In this work we want to investigate the importance of the DC mechanism for 
(n, 7 )-reactions on intermediate- and heavy-mass target nuclei occuring in the s- 
and r-process. In Section 2 we introduce the DC formalism and the folding proce- 
dure. In Section 3 the relativistic mean field theory (RMFT) and its application to 
neutron rich nuclei will be discussed. Finally in Section 4 we compare the influence 
of DC to radiative capture cross sections in the s- and r-process. 



2 Direct capture and folding procedure 

The theoretical cross section a th is obtained from the DC cross section cr DC given 

by 1, §, § 

o^ = Y.C?S t o? c . (1) 

i 

The sum extends over all possible final states (ground state and excited states) in 
the final nuclei. The isospin Clebsch-Gordan coefficients and spectroscopic factors 
are C{ and Si, repectively. The DC cross sections af c are essentially determined by 
the overlap of the scattering wave function in the entrance channel, the bound-state 
wave function in the exit channel and the multipole transition-operator. For the 
computation of the DC cross section we used the direct-capture code TEDCA ||. 

The folding procedure is used for calculating the nucleon-nucleus potentials in 
order to describe the elastic scattering data and the bound states. This method 
was already applied successfully in describing many nucleon-nucleus systems. In 
the folding approach the nuclear density p A is derived from experimental charge 
distributions |11] and folded with an energy and density dependent NN interaction 

V(R) = XV F (R) =\J p A (f)v cS (E,p A , \R - f\)dr . (2) 

with R being the separation of the centers of mass of the two colliding nuclei. The 
normalization factor A is adjusted to elastic scattering data and to bound- and 
resonant-state energies. The potential obtained in this way ensures the correct 
behavior of the wave functions in the nuclear exterior. At the low energies considered 
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in nucleosynthesis the imaginary parts of the optical potentials are small. The folding 
potentials of Eq. ^ were determined with the help of the computer code DFOLD 1C ] . 



3 Application of RMFT to neutron rich nuclei 

The RMFT describes the nucleus as a system of Dirac nucleons interacting via 
various meson fields. In the last few years this theory has turned out to be a very 
successful tool for the description of many nuclear properties (for example binding 
energies and charge radii for stable isotopes) p2| . The RMFT is built upon two 
main approximations: the mean-field approximation and the no-sea approximation. 
The mean-field approximation removes all quantum fluctuations of the meson fields 
and uses their expectation values. This approach cuts down all many body effects 
because the nucleons move as independent particles in the meson fields. The no-sea 
approach neglects all contributions from antiparticles. Only negative energy states 
are taken into account. One great advantage of the RMFT is that the spin-orbit 
interaction is described in a proper way without any additional parameters. 

The theory starts with an effective interaction of Dirac nucleons with mesons 



and the electromagnetic field. We work with the (au) /?)-model 12]. The cr-meson 
mediates the medium range attraction between the nucleons. The isoscalar vector 
mesons to cause a short-range repulsion. The contribution of the p-mesons is impor- 
tant for neutron- and proton-rich nuclei. There are six parameters which are usually 
obtained by fits to finite nuclear properties: each coupling constant of the meson 
fields with the nucleons g a , g w and g p , the constants of the nonlinear cr-potential (72 > 
53, and the mass of the cr-meson. The other meson masses are given empirically. 

In our calculations we have used the parameter sets NLSH (13], Jl5| and NL1 |lq| . 
The parameter set NL1 is fitted to the ground state binding energies and charge radii 
of spherical nuclei. This set overpredicts the neutron radii of neutron rich nuclei. 
The set NLSH has been fitted not only to binding energies and charge radii, but also 
to the neutron radii of several spherical nuclei. The binding energies of Sn-isotopes 
as a function of mass number show a kink at A = 132 (upper part of Fig. 1) signifying 
the shell effect at the magic neutron number TV = 82. The neutron-skin thickness of 
the Sn-isotopes also show this kink by about 0.25(A r — Z)/A (lower part of Fig. 2). 
The values for the neutron-skin thickness calculated with the parameter set NLSH 
are about 20% less than the values determined with NL1. The results of NLSH for 
the binding energies and neutron-skin thickness for stable neutron rich nuclei are 
comparable to the experimental data (see Fig. 1). For nuclei in the r-process (e.g. 

132 Sn _138 Sn) 

the neutron skin has a value of about 0.5 fm. 
With the RMFT we also can determine the proton and neutron density distri- 
bution for the Sn-isotopes. We found acceptable agreement between the calculated 
proton-density distributions obtained from RMFT for the stable Sn-isotopes 112 Sn- 
124 Sn with the experimental data [11]. 
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Figure 1: Upper part: Binding energies of Sn-isotopes calculated with RMFT using 
the parameter set NLSH. Lower part: Neutron-skin thickness of Sn-isotopes calcu- 
lated with RMFT using the parameter sets NL1 (broken curve) and NLSH (solid 
curve) compared with experimental data [17]. 
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Table 1. DC-calculations of radiative capture cross sections on neutron-magic tar- 
get nuclei compared with the experimental data. 



Reaction 



48 Ca(n, 7 ) 4y Ca 
86 Kr(n, 7 ) 87 Kr 
88 Sr(n, 7 ) 89 Sr 
136 Xe(n,7) 137 Xe 
138 Ba(n,7) 139 Ba 
208 Pb(n, 7 ) 209 Pb 



(a DC )(25keV) [mb] 



0.96 
0.09 
0.26 
0.16 
0.44 
0.13 



(a E xp)(25keV) [mb] 



1.03 ±0.09 
3.54 ±0.25 



23] 
Ml 



7.0 [25| 



1.05 ±0.09 |4| 
4.46 ±0.21 |26| 
0.31 M 



4 Direct capture for intermediate— and heavy— mass nu- 
clei 

In this section we want to investigate the importance of DC for intermediate- and 
heavy-mass target nuclei in the s- and r-process. The DC will be large compared 
with the CN contribution, if the level density of the compound nucleus is low, because 
then there will be only a few states in the compound nucleus, which can be excited 
in the reaction. This is true for example for the radiative capture reactions on light 
nuclei cited in the introduction. For intermediate- and heavy-mass nuclei this is 
also the case if the Q-value of the reaction as well as the projectile energy are small. 

For radiative capture on intermediate- and heavy-mass target nuclei induced 
by charged particles in astrophysical scenarios the projectile energy is so large that 
DC can be neglected. For instance, the DC-contribution of 144 Sm(a,7) 148 Gd is at 
least 5 orders of magnitudes lower than the CN-contribution However, a low 
level density of the compound nucleus can occur for radiative capture of neutrons 
on neutron-magic and/or neutron rich target nuclei in the s- or r-process. 

In Table 1 we show the results of the DC-contributions of (n,7)-reactions for 
some neutron-magic target nuclei with N = 28 (2(J, N = 50,82 [gl| and N = 126 
|22^ 1 occuring in the s-process. The DC-calulations were performed using the folding 
procedure for the optical and bound-state potentials as presented in Sect. 2. The 
average cross sections < a > at kT = 25 keV are compared to the experimental data. 
As can be seen the DC gives non-negligible contributions to the total cross sections. 
A special case is neutron capture on the double magic, neutron rich nucleus 48 Ca, 
where up to about 1 MeV no compound-nucleus levels exist, which can be excited. 
Therefore, for this reaction the cross section is almost given entirely through DC. 

As an example for neutron rich nuclei we investigated the radiative capture on 
Sn-isotopes up to the r-process path. The cross sections for such (n, 7 )-reactions are 
needed for the description of the freeze-out in the r-process [^7| . Since in this case 
no experimental data is available we compared the DC- with the CN-calculations 
as obtained from the statistical Hauser-Feshbach (HF) method. In order to make 
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Table 2. DC- and HF-calculations of radiative capture cross sections on neutron 
rich nuclei. 



Reaction 


Q-value [MeV] 


(cj D c)(30keV) \jjb] 


(a HF )(30keV) [/ib] 


i ^Sn(n, 7 ) i33 Sn 


2.581 


220.7 


116.4 


134 Sn(n, 7 ) 135 Sn 


1.871 


132.4 


32.5 


136 Sn(n, 7 ) 137 Sn 


1.611 


101.0 


18.5 


138 Sn(n, 7 ) 139 Sn 


1.230 


60.4 


8.4 



a meaningful comparison between the DC- and CN-reaction mechanism, we used 
the same masses, Q-values, spin assignments and excitation energies and optical 
potentials in both calculations. 

The densities necessary for the determination of neutron-nucleus folding poten- 
tials (Eq. 2) involving unstable nuclei cannot be taken any more from experimental 
data. We obtained them from the RMFT using the parameter set NLSH as described 
in Section 3. The strengths of the folding potentials were adjusted to reproduce the 
same value of the volume integral of 425 MeV fm 3 as determined from the experi- 
mental elastic scattering data on the stable Sn-isotopes |2lj[] , p9p. In order to obtain 
experimentally unknown masses and Q-values a microscopic-macroscopic mass for- 
mula based on the FRDM was used p0| . The spin assigments and the excitation 
energies of the ground and low-excited states of the odd residual nuclei 

133 Sn „139 Sn 

have been taken from [31]. For the neutron spectroscopic factors necessary fo the 
DC-calculation a value of unity was assumed. This assumption should be reasonable 
for the neutron rich Sn-isotopes involved in the r-process. 

The statistical-model calculations were performed with the code SMOKER 32]. 
Above the highest known state a level density description based on the backshifted 
Fermi gas |33| is employed, with parameters as given in |34[. The imaginary parts 



of the optical potential necessary for the HF-calculations were taken from [35]. 

In Table 2 we show the results of the DC- compared to the HF-calculations for 
the averaged cross section at 30keV on the even target nuclei 132 Sn- 138 Sn. As can 
be seen from this table both cross sections decrease when going to lower Q-values. 
However, the HF-cross section decreases much faster due to the lower level density 
of the compound nucleus. For the neutron-capture on 138 Sn the DC is about 7 times 
as high as the CN contribution as obtained with the HF-method. 

Summarizing, we showed that DC is not only important for radiative capture on 
light nuclei, but is also non-negligible for magic and neutron rich target nuclei. For 
some double magic and neutron rich nuclei in the r-process DC is even the dominant 
reaction mechanism. 
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